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ABSTRACT 

A dispersion relation involving all space -time harmonics has 
been obtained for electromagnetic wave propagation in an infinite 
dielectric njiedium whose permittivity is modulated by a high inten- 
sity pump wave containing two harmonics. Brillouin diagrams are 
given for the parametric coupling between the n = -1 (difference 
frequency) and n = +1 (sum frequency) space -time harmonics of the 
electromagnetic signal wave. The effect of the amplitude and phase 
of the second harmonic of the pump wave on the parametric growth 
process has been studied. Both stable and temporally unstable inter- 
actions are considered. 



Traveling wave parametric interactions of electromagnetic 
waves propagating in a dielectric medium with permittivity varying 
periodically!- in time and space has assumed considerable importance 
because of its numerous practical applications. Although many 
theoretical papers [l]-[4] have appeared on this problem, it is as- 
sumed for the most part that the perturbing pump wave modulating 
the medium is a weak, progressive, sinusoidal disturbance. However, 
there are a number of physical situations in which the modulating pump 
wave can have a large amplitude and contain two or more harmonics 
of the fundamental pumping frequency. Stimulated Brillouin scattering 
processes in Quantum Electronics [5], optical diffraction by non- 
linear ultrasonic waves [6] or electromagnetic diffraction by large 
amplitude traveling wave disturbances in the ionosphere or in labora- 
tory plasmas [7] are all various manifestations of this phenomena. 

The purpose of this report is to study the effect of the amplitude and 
phase of the second harmonic component of the pump wave on the para- 
metric interactions resulting from electromagnetic wave propagation 
in dielectric media of this type. 

It is assumed that the pump wave traveling along the z direction 
in an infinite, dielectric medium contains two harmonics and causes 
the dielectric constant of the medium to vary in the following manner 

e(z, t) = e Q + cos(w v t - k v z) + e 2 cos(2w v t - 2k y z + *0 (1) 

where e 2 « Cj « €q. W v and k v are the angular frequency and wave 
number of the fundamental component of the pump wave and fi is the 
phase angle of the second harmonic relative to the fundamental. The 
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pump wave at frequency 2u v can either be the second harmonic of a 
high-intensity, non-linear pump wave in which case the phase angle 

j 

= 0, or an! independent wave generated by a separate source in 
which case the phase angle 4 can be retarded by various values. For 
the purpose of this report only the phase angles ^.= 0 and = it are 


considered here. 


In Equation 1 it is assumed that 

€■» j« e, « e A . 


This assumption is true only if the punrip wave has a weak second har- 
monic component; in most of the physical situations described earlier 
the second harmonic pump component is indeed small and Equation 2 
is applicable [5] -[7]. 

Let us assume a TEfoi type electromagnetic wave moving in the 


z direction and parallel to the pump wave. The wave equation for the 


electric field E(z, t) in the medium can be written as i 


its 2 8t 2 *’ 

Pq is the peimeability of the medium and the displacement D .= e(z, t)E 
where e(z,t) is given by (1). By using Floquet’s theorem in conjunction 

with the principle of superposition, it can be shown that the electric 

! 

field can be represented in the following form: 


E(z,t) = E n e 




-jn(w t-k z) 

J ' V V 


where and k^ are the angular frequency and wave number of the signal 
wave in the unmodulated medium. The phase angle = 0 will be considered 
first and the dispersion relation for the $ - It case will be obtained by 
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making the necessary change in the sign of the second harmonic 
modulating coefficient later in this report. By substituting (4) 
in (3) we obtain the following fourth-order, linear, homogeneous 
difference-differential equation; when the phase angle f& = 0: 


D A + [A 
n n 2e Q L 1 


T + A ,, 
n-1 n+1 


Jf2_ 


] + x* t A n -2 + A n+2^ = 0 


(5) 


where n - 0, tl, £ 2, . . . and 


D = 
n 


1 - c 


2 ( k l + ^y) 

0 (u £ + nw v ) 2 


(6) 


c 0 = 1 /^o e o is the velocity of the electromagnetic wave in the unmodu- 
lated medium. Eq. (5) represents an infinite set of linear difference 
equations in the modal amplitudes A . It will now be shown that solutions 
can be obtained in the form of rapidly converging continued fractions. 

By adding the (n-l)th and (n+l)th equations given by (5) and sub- 
stituting the resulting expression for [A n _^ + A n ^.2^ * n we obtain: 

e, e 2 e ? 

DA + ~^[A .+ A. ,]--=■ [D .A .+ D ,,A A 

n n 2 cq l n-1 n+l J n-1 n-1 n+1 n+l J 6 q n 



lA 


■ i + A , t 
n- 1 n+1 


M 



- + A , _] 
n-3 n+3 J 


= 0 (7) 

From Eq. (2) it is seen that (c^/cq) <-< 1 and « 1. Hence, the 

coefficients of the last two terms in (7) which involve the quadratic 
products of (e^/^o) and ^2.^]) are mu °b smaller than the coefficients 
of the other terms in (7). Therefore, the last two terms in (7) can be 
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neglected and we obtain a linear, second-order difference-differential 
equation with rational coefficients: 



» -,]A , +[' 


The electric field amplitudes can be expressed now in the form of con- 
tinued fractions as follows: 



and 



Equations 9 and 10 apply to the case where 4 r The electric field 
amplitudes when the phase angle 4 - If can be obtained by replacing 
in the equations 5 to 10 by -e£. 

til 

A dispersion relation for the n space-time harmonic in the form 
of continued fractions, for the two cases when 4 ~ 0 and 4 - 7T is given by 
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The te^ in the above equation indicates whether the second harmonic 
of the pump wave is in phase ( </> = 0) or out of phase = Tf) with the 
fundamental component of the pump wave, respectively. The conver- 
gence of the continued fractions in (9), (10) and (11) will now be proved. 
From (8) and (6) it is seen that the following limit exists: 


Lt D 

n 

n-*co 


Lt D . , = Lt D , 
n+1 n-1 

n-*oo n-*oo 



( 12 ) 


where V = w v /^- v = the phase velocity of the pump wave. 

In the asymptotic limit n— >oo the coefficients of the difference- 
differential equation (8) are constants. Hence, by Laplace's method, it 
can be now shown that equation (8) is a Poincare difference equation [3, 8] 
whose characteristic equation is given by: 
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By applying Poincare's theorem [3, 8] it can be shown that the finite limits 

A - * A * 

Lt -gii = Pj and Lt = P 2 exist if 

n~»oo n i n-»oo A n 



(14) 


When (e^/ej) 
limits except 


« 1, the above condition indicates the existence of finite 
in the sonic region given by 


/i 

V ! e o 


^ V ^ 


where v = VS 


[1 +”] 
e l 


f?. 


(15) 


[1 ~J 


2/2 


/cq. The existence of a solution for (8) automatically 


assures the convergence of the Continued fractions in (9), (10) and (11). 


If the second 
then eqs. (8) 


harmonic modulating coefficient e 2 is set equal to zero 
■(11) and (14) reduce to the form given by Hessel and 
Oliner [9}. In particular, when e 2 = 0, (15) reduces to the form 

„ _i .. < 1 

fW. FZ 


This conditiok is identical to the limits of the sonic region defined by 
Cassedy and Oliner [1] for wave propagation in a sinusoidally modulated 
dielectric medium. 

Outside the sonic region defined in Equation 1 5, the parametric 
interactions can be divided into two types as follows: Type 1 . When 
the phase velocity of the pump wave is less than the phase velocity of 
the signal wave (v = V/cq < 1), then the parametric interactions are of 
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the stable, non-inverting modulator type with the dispersion diagram 
showing stop bands in the phase synchronous region. Stimulated 
Brillouin scattering and optical diffraction by ultrasound come under 
this category^ Figs, la and b show the Brillouin diagram (u vs. k) 
for the coupling between n = -1 (difference frequency) and n = +1 (sum 
frequency) space -time harmonics . Stokes and Antistokes interactions 


of this type are of interest in non-linear optics and in ultrasonics which 
was one of thb main reasons for considering them in this report. The 
parameters chosen in the calculations were = 0. 5, v = V/c^ = 0. 25, 

and 0. 1, and -0. 1. The solutions were obtained with an IBM 

7094 using a jrap'idly iterative secant technique for solving the trans- 
cendental equation. It can be seen from Fig. 1 that the stop band region 
expands and Contracts depending on the amplitude and phase of the second 
harmonic purkp modulating coefficient In the absence of the second 

pump harmoriic the stop band region in terms of k^X^ extends from 
1. 974# to 2. 0987T. When the second harmonic is out of phase with the 
fundamental component of the pump wave = ir) and | e^/e^ | = +0. 1, 
the width of the stop band increases and extends from 1. 9531T to 2. 1237T. 
When the second harmonic is in phase with the fundamental = 0) and 
| j = 0. 1, the stop band width decreases to 1. 9 9 3 7T to 2. 077ir. 

Within the stop band, the solutions to the dispersion relation re- 
quire compleix propagation constants indicating contraflow directional 


coupling resulting in frequency down conversion when the signal and 
pump waves are traveling in the same direction. For each set of values 
of e£, and j i the real part of the propagation constant (3^ i® constant 
throughout the stop band while the imaginary component varies in a 
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semicircular manner. The complex propagation constant within the 
stop band is also significantly affected by the second harmonic pump 


wave. When = the maximum value of the normalized complex pro- 
pagation constant within the stop band is given by 0. 5317T + jO. 06 77T$ 
with the second harmonic out of phase = ir) and | = + 0 . 1 , 

the maximum value of the propagation constant in the stop band increases 
to 0. 5397T + jO. 0931T and reduces to 0. 523JT + jO. 0447T for the in-phase 
condition { ft = 0, I e 2.^ e l J ~ *)• From these computations it appears 

that the parametric coupling between the n * +1 and n = -1 space -time 
harmonics of the signal wave is stronger when the second harmonic of 
the pump wave is driven out of phase with the fundamental and decreases 
when it is in phase with the fundamental component of the pump wave. 


Type 2i . When the phase velocity of the pump wave is larger than 
the phase velocity of the signal wave in the unperturbed medium (v > 1 ), 
the parametric interactions are characterized by 'temporal instabilities' 


resulting in time growing oscillations. The nature of this oscillatory 
growth process has been well described by Cassedy [ 2 ]. The stop bands 

i - . 

in the Brillouin diagram are inverted, wherein solutions for the disper- 
sion relation are obtained only for complex frequencies and real wave 
numbers. The amplitude and phase of the second harmonic once again 
has a strong j influence on both the width of the stop band in temporal 
frequency (k 0 = w/c) and on the growth rate of the oscillations; however, 
the effect of the phase relation of the second harmonic on the parametric 


growth rate is opposite to that observed in the Type 1 interaction where 
v < 1 . Figs. 2a and b show the Brillouin diagram for the interaction 
between the n = -1 and n = +1 harmonics. When = 0 , the inverted stop 
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band extends from P T X = 4. 6l0?r to 4. 720 IT. When - 0 (in phase) 
and = 0. 1, the stop band expands in width and covers the range 

= 4. 625f to 4. 785 jt. For fixed values of ej, and the 
temporal frequency within the stop band is constant but the growth rate 
(imaginary part of k^X^) varies in an approximately semicircular manner. 
The effect of the second harmonic pump wave on the growth rate is shown 
in Fig. 2b* when = 0, the maximum value of the complex k^X^ in the 
stop band is 2. 25071 + jO. 042fl\ When | e^/e^ | = +0. 1 and $ = 0 (in phase 
condition), the maximum value of k^X^ is 2. 232?r + jO. 09^* The growth 
rate of the instability, therefore, increases with the amplitude of the 
second harmonic when it is in phase with the fundamental component of 
the pump wave. 
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